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ABSTRACT 

Observations of line emission from water masers near the center of the galaxy 
NGC 4258 have recently provided compelling evidence for rotating disk of gas, 
viewed nearly edge-on, surrounding a massive black hole. We show that the disk is 
very likely to be only marginally stable to radial perturbations - a stability regime 
where weak non-axisymmetric disturbances grow via the "swing amplification" 
effect, leading to the formation of a ragged, multi-armed spiral pattern similar to 
that observed in Sc galaxies. This suggests a natural explanation for the apparent 
clustering of the high- velocity emission sources into several distinct clumps, and 
for the observed regularity in the distance intervals between them. The clumps 
of maser sources appear at the intersections of the spiral arms and the radial line 
of longest coherent gain path (the diameter through the disk perpendicular to 
the line- of- sight), and are thus spaced apart at the characteristic crest-to-crest 
radial distance between the arms. 

This interpretation implies a disk thickness of ~ 0.003 (Q/l-G) pc at a radius 
of ~ 0.2 pc, where the local value of the stability parameter is 1.2 ^ Q ^ 2. 
The H2 density is ~ 1.8 X 10^° (Q/1.6)~^ cm~"^, assuming that H2 dominates the 
mass density in that region. The disk mass is ^ lO^M©, which is consistent 
with the accuracy of the Keplerian fit to the rotation curve of the maser 
emission sources. These results closely concur with previous estimates based 
on independent considerations. 

Subject headings: masers - accretion disks - galaxies: individual (NGC 4258) - 
galaxies: nuclei 
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1. INTRODUCTION 



The maser line emission from the circumnuclear disk in NGC 4258 consists of two 
groups of sources (Miyoshi et al. 1995). The systemic features appear along the line-of- 
sight to the center of the disk and have velocities comparable to the systemic velocity of the 
galaxy. The high- velocity features are distributed in a nearly planar structure at distances 
between 0.16 pc (hereafter r^) and 0.25 pc {rout) from the center of rotation, and their line- 
of-sight velocities are described by Keplerian motion to high precision. This remarkable 
finding confirmed earlier suggestions (Watson & Wallin 1994) that the high- velocity sources 
lie along the diameter through the disk (mid-line), perpendicular to the line- of- sight. The 
restricted spatial distribution of masers along the line-of-sight to the nucleus and along the 
mid-line of the disk are undoubtedly related to the fact that in these regions the gradient 
of the line-of-sight velocity is zero and hence, these regions support long path-lenghts for 
coherent amplification in our direction (Miyoshi et al. 1995). Obviously, all observers in the 
universe who view the disk nearly edge-on would see a qualitatively similar picture. 

The high- velocity maser sources are not distributed uniformly nor randomly along the 
mid-line. They are clustered into seven distinct "clumps," five on one side of the disk, and two 
on the other (figure 2a in Miyoshi et al. [1995]). Moreover, there seem to be some regularity 
in the distance intervals between the clumps; they are almost equally spaced apart from 
each other, with mean inter-clump separation of ~ 0.023 pc (0.75 mas). It is plausible that 
the intriguing clustering of sources is due to small gas clumps in the disk inside which the 
conditions for amplification of maser emission are favorable for as yet unknown reason, but 
this simple interpretation is improbable for the following reason: if noticeable maser emission 
had been produced due to amplification over the column density within a gas clump, the 
chances of detecting that emission would have been independent of the location of the gas 
clump in the disk. In such case, we should have detected emission from gas clumps also at 
locations other than along the mid-line - a situation which is firmly ruled out by the observed 
nearly perfect Keplerian rotation curve. 

We shall now argue that the observed clustering of the maser sources, and the indicated 
regularity in the distance intervals between them can be naturally explained by simple disk 
stability considerations (§2). In §3 we conclude and discuss caveats. 



2. DISK STABILITY AND SPIRAL STRUCTURE 



2.1. The Condition For Spiral Structure 



A rotating, self-gravitating cold disk is locally stable against axisymmetric disturbances 
of wavelength A > Xcrit^ where 

_ 47r2GS 

^crit — 
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S is the surface density of mass per unit area, K,(r ) = 0[4 + c/log(02)/(/logr]''' is the 
local epicycle frequency, and 0(r) is the angular velocity. A non-rotating sheet with non- 
zero sound speed, u^, is stable against perturbations with A < u^/GS. Thus, under the 
joint influence of rotation and moderate superposed sound speed, the unstable axisymmetric 
disturbances are confined to a certain intermediate range of scales - a range that shrinks 
as the sound speed is increased. A thin gaseous disk is locally stable to all axisymmetric 
perturbations when Q > 1^ where (Toomre 1964; Julian & Toomre 1966) 



Although Q > 1 curbs the axisymmetric Jeans instabilities, a considerable growth of 
non-axisymmetric disturbances occurs when 1<Q^2^ leading to the formation of a multi- 
armed spiral density pattern (Goldreich & Lynden-Bell 1964; Julian & Toomre 1966; Toomre 
1981, 1990). This phenomenon was discovered in the context of spiral galaxies (Goldreich 
& Lynden-Bell; Julian & Toomre 1966), and is believed to be the basic engine behind the 
structure of Sc galaxies. In contrast to "grand design" spirals like M51 or M81, Sc galaxies 
exhibit a ragged pattern of many spiral arms. The underlying physical mechanism which 
is responsible for this structure is the "swing amplification." It results from a three-fold 
conspiracy between the shear flow, epicyclic shaking and self-gravity (Goldreich & Lynden- 
Bell 1964; Julian & Toomre 1966). The disk particles travel almost freely through the density 
pattern, so any individual spiral arm does not remain composed always of the same particles. 
The spiral pattern is constantly drifting and shearing with the differential rotation, and spiral 
arms are constantly forming and dying with typical life-time of order half the revolution 
period. As observed in simulations (Miller, Prendergast & Quirk 1970; Hohl 1970; Hockney 
& Brownrigg 1974; Toomre 1981; Sellwood & Carlberg 1984; Carlberg & Freedman 1985), 
even weak transient perturbations with relative force of a fraction of a percent can yield 
an evolving spiral pattern of impressive severity. The maximum growth occurs for waves of 
statistically well-defined wavelength (Goldreich & Lynden-Bell 1964; Julian & Toomre 1966), 
leading to a spiral structure which is most pronounced at radial crest-to-crest separation of 
~ Acrit, as defined in equation (1). The less self-gravitating the disk, the smaller is the radial 
distance between the arms, and the larger is the number of arms {e.g.^ Carlberg & Freedman 
1985). 



2.2. The Circumnuclear Disk 

We suggest that the stability state of the circumnuclear disk in NGC 4258 is similar 
to that of Sc galaxies, implying that it too may well have a ragged, multi-armed spiral 
structure within some range of radii. This could explain the apparent clustering of maser 
emission sources as appearing where spiral arms intersect the mid-line (Figure 1). Since 
these intersections are spaced ~ Acrit apart, it would produce some regularity in the distance 
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intervals between the clumps, as observed. We shall now derive the disk characteristics 
required for the existence of a spiral structure at the region of the high- velocity sources. 

The mean distance of the high- velocity sources from the center of rotation is f ~0.20pc. 
The thickness of the disk, ^, at radius f, and the local circular velocity, are related by the 
condition of an hydrostatic equilibrium which reads ^/f ~ Ug/u, where Vg is the effective sound 
speed in the disk at radius f (notice that Vg need not necessarily be directly related to the 
temperature of the masering gas, e.g.^ in the case of the interstellar medium in galactic disks 
the internal temperature of diffuse Hi clouds is significantly lower than the root mean-square 
random cloud velocity [Spitzer 1978]). Since the force field is Keplerian to high precision 
(Miyoshi et al. 1995) we have n{f) = 0(f) = vjf = 1.42 X 10~^° s~^. Combining the above 
equalities with equations (1) and (2) we obtain 

h = 9^ ~ 0.1 (0/1.6) mas , (3) 

where 1.2 ^2 is required for the existence of a spiral structure, and Xcrit — 0.75 mas 
is the mean radial distance between the clumps. Assuming that H2 makes the dominant 
contribution to the mass density of the disk at radius f, and that the relative fraction of 
water vapor is negligible {e.g.^ Anderson & Watson 1993), the density of H2 molecules at f, 
averaged over the thickness of the disk, is 

= ^, ^""^/^ - 1.8 X 10^° (0/1.6)"' cm-3 . (4) 

For comparison, Miyoshi tt al. (1995) estimated a disk thickness of < 0.1 mas and H2 density 
of ~10^°cm~"^ based on independent considerations. 

Since the disk is not massless it must introduce some systematic deviation from a perfect 
Keplerian rotation curve. Although the surface mass density profile of the disk is unknown, 
the above derivation fixed its value at radius ~f [J] = nhniYi.^) . Let us assume, for example, 
that S(r) = J]f/r. The rotation profile in the absence of the black hole would have been flat, 
with u^ = 27rGSf. Thus, the rotation curve in the combined force field is v^(r) = (a/r) + u^, 
rather than v'^{r) = a' jr. The systematic deviation of v{r) from v'{r) over the entire range 
of radii [ri„, rout] is given by ~ v(ri„) — v'(ri„) under the constraint that v(rout) = v'{rout)- 
We obtain 

. 1.8 X 1010 cm-3y l^o.lmas^ 

which is independent of Q. Such deviation is consistent with the accuracy of the 
Keplerian fit to the rotation curve (Jim Herrnstein, private communication). The minus 
sign implies that the rotation curve should be slightly flatter than an r~^^^ power law. 
If S(r) oc all the way to the center, the total mass of the disk within Vout is 

~ 9 X 10^ (n/1.8 X 10^° cm~"^) (^/O.l mas) M©. Assuming a steeper surface density profile 
at f, in which case the associated rotation curve falls with distance, would result in even 
smaller deviation from a Keplerian relation. 



A- - -4.2 ( ...,:S.^-3 )( 777^1 kms- , (5) 
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3. DISCUSSION 

We have shown that the puzzHng distribution of high- velocity maser emission sources 
can be attributed to a spiral structure in the circumnuclear disk. This interpretation implies 
a mean disk thickness (Eq. [3]) and density (Eq. [4]) which closely concur with previous 
estimates that are based on independent considerations. While spiral structure on scales of 
~ 10'* parsec in Sc galaxies is outlined by tracers such as Hn-regions, OB-star associations, 
21-cm line emission, and CO lines, a similar spiral pattern on a sub-parsec scale appears to 
be outlined by maser line emission, thus providing unprecedented insight into the structure 
of circumnuclear disks. 

The spiral structure interpretation implies that the conditions for amplification of maser 
emission are more favorable in the wave crests than in the troughs. This raises an interesting 
question: how large should the density contrast in the spiral pattern be, relative to the mean 
density of the disk, in order for the maser emission to be produced preferentially in the 
spirals arms? 

Answering this question is not straightforward since the condition for the onset of 
masering, as well as the maser luminosity, depend not only on the density but also on 
the temperature, the H2O /H2 ratio, the amount of dust grains {e.g.^ CoUison & Watson 
1995), the nature of the pumping mechanism, and on whether the emission is saturated or 
not. There is currently no theory which gives the maser luminosity as a function of all these 
factors (M. Elitzur, C. McKee & D. HoUenbach, private communication), except in a few 
specific cases. For example, assuming that the pumping mechanism is coUisional excitation, 
which is likely to be the case for the high- velocity masers, the brightness temperature of 
the maser emission has been shown to scale as /Og^^, where pgas is the gas density (Elitzur, 
HoUenbach & McKee 1989, 1995 in preparation). Given this non-linear dependence on 
density, and the fact that the S/N ratio for the intensity of the high- velocity features is ^10 
(Figure 1 of Miyoshi et al. 1995), it is clear that a density ratio of p crest I P trough ^3 would 
be enough to reproduce the observations (notice that p crest I P disk < P crest / P trough)- A density 
contrast of order unity can be easily attained by a spiral pattern. It has been shown (Toomre 
1981, and references therein) that the density contrast of a spiral arm depends on the local 
value of the Q parameter, and it changes in time as the arms shear, wind up, and disperse. 
The maximum amplification factor varies between ~3 when Q — 2, to ~ 100 when Q — 

In fact, the observed distribution of maser sources can be reproduced even by a lower 
density contrast than that suggested in the above discussion. The gas temperature in the 
spiral arms is likely to be higher than that between the arms, and it is known that the 
brightness temperature increases (in a non-trivial way) with temperature (D. HoUenbach, 
private communication). Studying the dependence of the maser luminosity on all the 
potentially important factors is far beyond the scope of this Letter and is left for a future 
investigation. Again, we emphasis that the condition for spiral activity and the radial 
distance between the arms are independent of the density contrast. 

The existence of a spiral structure in the disk implies that the velocity field must deviate 
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to some degree from a smooth Keplerian rotation curve. The ampHtude of the deviations 
from circular velocities (hereafter, the peculiar velocity) is comparable to the amplitude of 
the epicyclic velocity around the guiding center, and it can be estimated in the following way: 
the orbital radius of the epicyclic motion is ~ Acrit/2, where Xcrit is the mean distance between 
the arms. In the force field of a point mass the epicyclic frequency is equal to the orbital 
frequency of the guiding center and we have Vepicydic — {v^crit 1'^^) — 45km s~^, where v is the 
orbital velocity at radius f. Although such velocities are larger than the measured accuracy 
of the Keplerian rotation curve (~4kms~^), we should not be able to detect these motions 
for a simple reason: at any point in the disk the local peculiar acceleration due to the spiral 
pattern is always orthogonal to the wave troughs or crests (Toomre 1981), and the peculiar 
velocities are predominantly along the same direction. Since the high- velocity maser emission 
is observed only along diameters through the disk which do not deviate from the mid-line 
by more than 6°~0.1rad (L. Greenhill, private communication), the peculiar velocities are 
almost perpendicular to the line-of-sight and are thus almost undetectable (the component 
of the peculiar velocity along the line-of-sight would not exceed ~0. 1x45 = 4. 5kms~^). Since 
the detected deviations from a Keplerian rotation curve are not much larger than the above 
value, we can also conclude that the pitch angle of the spiral arms, within the observed range 
of radii, does not exceed ~6°. 

Explaining the difference in the distribution and intensity of red-shifted versus blue- 
shifted sources remains a serious challenge for theory. It is currently hard to understand this 
asymmetry, regardless of whether there is a spiral structure or not. It might be related to 
the fact that the disk is slightly warped (see Neufeld & Meloney 1995; Miyoshi et al. 1995). 

We should emphasize that finding the disk near the border of stability [1 < Q -^2) 
does not necessarily requires fine tuning (Goldreich & Lynden-Bell 1964; Julian & Toomre 
1966). Although the individual disk particles travel almost freely through the wave pattern, 
the constantly-evolving spirals tend to slowly increase the gas temperature, thus increasing 
the Q-parameter. In order to maintain a spiral activity the disk must be cooled. In Sc 
galaxies, energy dissipation in small regions in the gas leads to the formation of condensations 
such as giant molecular clouds, a process which effectively cools down the disk. The non- 
axisymmetric disturbances associated with the new condensations trigger new modes, and 
generation after generation of spiral arms form, wind up, and disperse. This has been 
suggested to be the self-regulating, regenerative mechanism of spiral structure in Sc's, and 
it may as well apply to the circumnuclear disk in NGC 4258. If the temperature of the disk 
is indeed maintained by this process, the implied slow depletion of gas should be balanced 
by continuous gas accretion. 

Finally, we point out an alternative conceivable explanation for the clustering of the 
high- velocity sources. It has been shown (Syer, Clarke & Rees 1991, and references therein) 
that the cumulative effect of many passages of a star through the disk would bring it into a 
circular orbit corotating with the disk. If the disk is thinner than the star's Roche radius, 
the star may open a gap in the disk. The clustering of the maser sources could then be a 
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manifestation of the ring-like structure of the disk due to the grooves produced by many such 
stars. However, we find this interpretation improbable mainly because the distance intervals 
between the clumps are substantially larger than the gravitational "sphere of influence" of 
a star in that disk. 

I thank Alar Toomre, George Field, Jim Moran, Lincoln Greenhill, Jim Herrnstein, 
Charles Gammie, Moshe Elitzur, and David HoUenbach for discussions. This work was 
curried out at the Harvard- Smithsonian Center for Astrophysics, and was supported by the 
U.S. National Science Foundation, grant PHY-91-06678. 
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Fig. 1. — An illustraition of a multi-araied spiral structure in the circumnuclear disk. The 
solid line represents the "mid- hue" along which the coherent gain path for ajuphRcation 
of maser emission in our direction is long. It is argued that the clumps of high-velocity 
sources appear at the intersections of the mid-line and the spiral arms (see text). It should 
be emphasized that the real spiral pattern is unliltely to be a£ regular a£ in the above 
illustration. It probably resembles the structure of Sc galajdes which can be described a& a 
swirling hotch-potch of pieces of spiral arms. The apparent deficiency in the number of blue 
shifted clumps may be attributed to the indicated large scale warp in the disk. 



